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We present a feasibility study for a new analysis for extracting the angle 7 of the Unitarity Triangle 
from the study of the neutral B meson decays. We reconstruct the decay channel B° — > D°K*° with 
the K*° — > K + tv~ and the D — » Ks^ + tv~ using a D° Dalitz analysis technique. The sensitivity 
to the angle 7 comes from the interference of the b — > c and b — > it processes contributing to the 
same final state and by the fact that the B°(B°) can be unambiguosly identified through the sign 
of electric charge of the kaon from K*°(K*°) decay. The impact of the result of such analysis is 
evaluated for the actual BaBar statistics. 



I. THE TECHNIQUE 

Various methods related to B~ — > D^°K^~ decays 
have been proposed to determine the UT angle 7. Within 
them, the one that gives the best error on 7 is the Dalitz 
method Q0. 

We propose to measure the angle 7 using a Dalitz anal- 
ysis applied to neutral B mesons. In general, since the 
neutral B mesons mix, interference effects between b — ► c 
and b —> u decay amplitudes in B° decays (for instance 
into D^'^ir^ final states) are studied for the determina- 
tion of the combination of UT angles 2/3 + 7. I n this case 
the tagging technique and a time dependent analysis are 
required. This can be avoided if the final states contain 
a particle which allows to unambiguously identify if a B° 
or B° has decayed. This is the case of neutral B mesons 
decaying into D°[K + ir~] final states, where the flavor of 
the B meson can be determined through the sign of the 
electric charge of the Kaon. 

One of the relevant parameter in those kind of analysis 
is the ratio between the b — > u and b — > c amplitudes, 
which can be expressed by : 



r B (D°K*°) = 



\A(D°K*°)\ = \V ub \ A stmng (D°K*°) 
\A{D°K*°)\ MVcb\A stmng (D°K*°)' 



The sensitivity to 7 is proportional to the tb value. 
Considering the CKM factors in the b — > c and b — > u 
transitions and the fact that in both cases the processes 
are mediated by a color suppressed diagram, we expect 
rB(D°K*°) to be in the range [0.3-0.5], larger than for 
the equivalent ratio in the charged B sector, which has 
been found to be of the order 0.1 ||. 



II. STRUCTURE OF THE ANALYSIS 

In the analysis on data we would perform an unbinned 
maximum likelihood fit in order to extract the interest- 



ing quantities: the likelihood will contain a product of a 
yields pdf and a CP pdf (only the latter has, within its 
variables, the weak phase 7). 

In our feasibility study, we assume the yields pdf to 
be a function of the classic BaBar analysis variable tues 
1 and of a variable that is able to discriminate between 
BB and continuum events (qq events, with q — u,d, s, c). 
For the distributions of those variables we make some 
realistic assumptions inspired by BaBar analysis. 

In writing the CP pdf, we have an additional difficulty, 
with respect to the B~ — » D°K~ Dalitz analysis. This 
is due to the fact that the natural width of the K*° is 
not small (~50 MeV) and hence the interference with 
the other B° — > -D (^7'')n On _K* processes may not be 
negligible. We follow the formalism introduced in 0] and 
also used in the B~ — > D°K*~ BaBar analysis in order 
to solve this problem. 

We introduce the effective quantities rg, k and 8s de- 
fined as : 



T(B° -> D°(K+n-)) _ Jdp\Al p \ 
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where < k < 1 for the Schwartz inequality and 8s S 
[0, 2_7r]. In the limit of a B ^2-body decay, such as B° — > 
DK° we have : rs — > rg, 83 — ► 8b (the phase difference 
between the V u b and the V c b amplitudes of EqQ} and 
k=l. A cp , A up are real and positive. The subscript c and 
u refer to the b — > c and b — > u transitions, respectively. 
The index p indicates the position in the phase space of 
D°K + tt-. 

Applying this formalism, we write the CP pdf as fol- 
lows: 



I/tI 2 



T(B°(B°) -» £>[-> K s Tr-7r + ]K-TT + (K + n-)) oc 
r 2 s I /± 1 2 + 2kr s I f T 1 1 f ± I cos( J s + S D (ml , m\ ) T 7) ,(4) 
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1 the beam-energy substituted mass m^s = y (v^/2) 2 ~ Ps 2 > 
where the asterisk denotes evaluation in the T(4s) CM frame 
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0.0507 


1~ 


1.0 


0.4 


1.412 


0.294 
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0.3 


0.12 


1.4324 


0.109 


2+ 


0.15 


0.06 


1.717 


0.322 


1" 


0.2 


0.08 



TABLE I: List of mass (M), widths (W) and quantum num- 
bers of the resonances considered in our model. The last two 
columns present the chosen values of the coefficients cij for 
the Vcb and V u b transitions respectively. Note that the phase 
Sj are not indicated and their choice is arbitrary. 




FIG. 1: Distribution ofrs and k in a region within ± 50 MeV 
the nominal mass of the K*° (892) resonance (mic-n in the 
range [0.7159-0.8951] GeV 2 ). This distribution has been ob- 
tained by randomly varying all the strong phases (between [0- 
2n]), the amplitudes (between [0.7-1.3] of their nominal value 
(except those relative to K* which has been fixed to 0.4 and for 
the V u b amplitude of the D Sj 2(2573) ± which has been varied 
between [0.-2.] of its nominal value) given in Table\Tj). 



where 5d(to^, m±) is the strong phase difference between 
/(mj.,ml) and /(m^,mj_) and r s , k and 6s are defined 
in Eqs. ([2J) and ([3]). We have simplified the notation using 
f± = /(m|,m|) and f T = /(m|,m|). 

The Dalitz structure of the decay D° — ► Ksirn is well 
known and it has already been used in Dalitz analysis in 
the charged B sector [l| Q. 

III. THE k AND rs PARAMETERS 

We performe a study to evaluate the possible varia- 
tion of rs and k on the B Dalitz plot. For this pur- 
pose we use a B Dalitz model as suggested by the recent 
measurements [H], [f|. The model assumed for the de- 
cay parametrizes the amplitude A at each point k of the 
Dalitz plot as a sum of two-body decay matrix elements 
and a non-resonant term according to the following ex- 
pression : 

A c k (u k) e l5 ^ -^a.e^B^'Kr.sj+a,^ (5) 
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We consider a region within ± 48 MeV from the nominal 
mass of the K*°(892) resonance and we obtain the dis- 
tribution of rs, k and k x rs by randomly varying all the 
strong phases between [0-27r] and the amplitudes up to 
±30% of their nominal value (shown in table [TJ) . Only for 
the D S! 2(2573) ± resonance, the variation is up to ±200% 
of its nominal value. 

The distribution we obtain for rs and k are shown in 
figure U 

The results is that, in K* Q (892) mass region, rs can 
vary between 0.3 and 0.45 depending upon the values of 
the phases and of the amplitudes contributing in the K*° 
region. In absence of pollution we would have expected 
i*s = r B = 0.4. The distribution of k is quite narrow, 
centered to 0.95 with an error of 0.03. 

For this reason the value of k can be fixed to 0.95 and 
the effect of its variation (within a reasonable interval) 
will be considered a systematic error. 



IV. CP FIT STUDIES 

We perform an intensive toy-MC study on the CP fit, 
assuming the actual BaBar statistics (350fb _1 ). We as- 
sume to have 35 signal events, about 2000 continuum 
background events and 250 BB background events. 

A. Polar coordinates 

We first perform toy-MC studies in polar coordinates. 
In figure [2] we observe the known linearity problem: due 
to the dependence of the likelihood on rs, we tend to get 
from the fit a value of rs higher than the true (generated) 
one and consequently to underestimate the error on 7. 



r B : fitted value VS generated one r B : fitted value VS generated one 




FIG. 2: Distribution of r£ Ji vs r^^ N from toy MC for dif- 
ferent generated values. The toy-MC have been generated as- 
suming 350fb _1 (left plot) and 10 times the statistics (right 
plot). 

It is evident that the linearity effect (still visible in 
the results of the tests at high statistics) is not the only 
problem that affects our measurement. In fact, as it can 
be seen from the left plot, also for high values of Tg 
the results for vs r^ EN do not tend to the curve 
rg lt = rg EN (the black curve in the plots). We thus 
conclude that we cannot fit in polar coordinates (with 
rs, 7 and 6s floating) because of the two effects: the 
linearity and the low statistics and that is the second one 
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that dominates (at least for r$ 
be for our measurement). 



0.4, as it is expected to 



B. Cartesian coordinates 

We then perform a toy-MC study in cartesian coordi- 
nates pj. The use of these coordinates normally solves 
the linearity problem ( [l[ , Q ) . 

From the toy-MC the four variables appear to be bi- 
ased and show a non-Gaussian behaviour. In table HH we 
summarize the results of toy-MC where we generated for 
the yields the values we expect on 350/6 -1 (left column) 
and the results of toy-MC for 10 times the now available 
BaBar statistics (right column). As it can be seen in the 
left column, all the variables are biased and the a of their 
pulls are not compatibles with 1. This effect disappears 
at high statistics. We conclude that, with the now avail- 
able statistics, we cannot perform the measurement in 
cartesian coordinates. 



ICte stat 



„PULL 

„PULL 
<J X _ 

PULL 
„PULL 



-0.52 ±0.05 
0.82 ±0.03 
-0.07 ±0.05 
0.78 ± 0.04 
-0.18 ±0.05 
0.79 ± 0.04 
0.40 ± 0.05 
0.79 ± 0.04 



-0.04 ± 0.05 
0.97 ± 0.04 
-0.02 ±0.05 
0.99 ± 0.04 
-0.05 ±0.06 
1.01 ±0.04 
-0.03 ±0.05 
1.03 ± 0.04 



TABLE II: Pull distributions for cartesian coordinates on 
500 toy-MC with a signal background configuration similar 
to the one we find on data (left column) and for 10 times the 
statistics (right column). 



C. Measurement strategy 

We conclude that, with the now available BaBar statis- 
tics, the measurement strategy would be extracting from 
the fit 7 as a scan with respect to rs (i-e. by performing 
a likelihood scan on 7 for each value of rs with 7, S and 
all the yields parameters floated in the fit). 



In this way we would extract from data the maximum 
possible information for the available statistics: an infor- 
mation on 7 and not on rs- The 7 distribution obtained 
in that way, already very interesting on its own, would 
be very precious when combined with an experimental 
input for tb- The coverage tests on toy-MC generated 
with different values of r$ show that we have no bias in 
this fit configuration. 

V. IMPACT OF THE MEASUREMENT 

In this section we show, on one chosen toy-MC, what 
would be the impact of a real mesurement. In figure [3] 
we show the output of our fit: a likelihood scan of 7 with 
respect to rs- 



Likelihood scan tor y as a tunction or r 




FIG. 3: Negative Log Likelihood scan for 7 for different values 
of tb- The distribution is obtained on a chosen toy-MC. 



This distribution, when combined with a fake r$ mea- 
surement (rs = 0.35 ± 0.15) gives 7 with an error of 46°. 
A measurement of rs with such an error is what we could 
expect to obtain, if rs — 0.35, on 500fb _1 data from an 
analysis of B° — > D°K*° with the D° in flavor modes. 

Clearly, this is just a example on a chosen toy-MC and 
it is not necessarily representative of what we will find 
on data. Indeed we found, on toy-MC studies, that in a 
small percentage of cases we can be much less sensitive 
to 7. That problem has been found to depend on the low 
number of signal events (such that in some cases those 
events happen to be in poorly sensitive regions of the 
D° Dalitz plot) and it will therefore disappear with the 
increasing of the statistics. Still, our test shows that we 
could have a result on 7 that is competitive with the 
charged B Dalitz analysis (the analysis made on BaBar 
data gives today an error on 7 of 42°). 
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